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We study one-component fermions in chain lattices with proximity-induced superconducting gap 
and interparticle short-range interaction, capable of hosting Majorana fermions. By systemati¬ 
cally tracking various physical quantities, we show that topological states and topological phase 
transitions in the system can be identified by multiple signatures in thermodynamic quantities and 
pair-condensate properties, in good agreement with the known signatures in the ground-state energy 
and entanglement spectrum. We find the disappearance of the topological phase in a largely attrac¬ 
tive regime, in which the system undergoes a first-order transition between two topologically trivial 
states. In addition, the stability of the signatures against finite size, disorder, and inhomogeneity is 
analyzed. Our results provide additional degrees of freedom for the characterization of topological 
states with interaction and for the experimental detection of emergent Majorana fermions. 
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I. INTRODUCTION 

Exploring topological states of matter has become a 
rapidly developing field in condensed matter physics [1- 
3]. One intriguing topological state that exhibits both 
fundamental interest and practical application is the 
emergent Majorana zero modes or Majorana fermions [4- 
8], which are their own antiparticles and possess zero 
energy, in superconducting materials. The pursuit of 
the Majorana modes started from the study of a one¬ 
dimensional (ID) p-wave superconducting chain by Ki- 
taev [9]. This milestone has triggered several alternative 
schemes for the realization of Majorana fermions, such as 
p-\-ip superconductors [10-14], spin-orbit coupled super¬ 
conducting nanowires [15-18], chains of magnetic atoms 
on superconducting substrates [19-22], superconducting 
surfaces of topological insulators [23-26], superfluid he¬ 
lium 3 [27-29], and ultracold atoms [30-35]. Recently, 
a semiconductor nanowire with intrinsic spin-orbit cou¬ 
pling, external magnetic held, and proximity-induced su¬ 
perconductivity has become one of the experimentally 
promising platforms to host Majorana fermions [36-41], 
which appear on the edges of the wire and lead to a tun¬ 
neling conductance peak at zero voltage. Such a zero-bias 
peak has been observed as (indirect) evidence for their 
existence. 

In addition to the transport properties [42-54] , search¬ 
ing different signatures for the topological phase is an 
ongoing task for the investigation of Majorana fermions. 
From the experimental point of view, it not only pro¬ 
vides more evidence for direct detection of Majorana 
fermions but also helps rule out different physical causes 
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that result in the same transport behavior [55-58] . From 
the theoretical point of view, a comparison between var¬ 
ious quantities, even those typically used to describe 
Ginzburg-Landau-type phase transitions, such as suscep¬ 
tibilities and superfluid order, can provide useful informa¬ 
tion for characterizing topological phases and topological 
phase transitions, especially in interacting systems [59- 
61]. Interaction effects, which are unavoidably present 
in reality, may alter physical features of the topological 
phase and even change the topology [59-71]. For exam¬ 
ple, time-reversal symmetric Kitaev chains in class BDI 
change the topological invariant from Z to Zg when the 
interaction is turned on [71]. In an interacting system, 
Majorana zero modes at edges become many-body Ma¬ 
jorana wave functions [72, 73] and the degenerate ground 
states with two different parities are connected by these 
many-body Majorana zero operators. Such many-body 
phenomena attract broad interest, from the fundamental 
understanding of its nature to applications on quantum 
computation [74-77]. The study of multiple signatures 
shall provide a convincing series of tests to characterize 
the interacting topological phase diagram. Recent works 
have analyzed individual quantities for separate models, 
such as entanglement spectrum in a spin-orbit coupled 
chain with interactions [62, 66], compressibility [78] and 
spectral function [63] in the Kitaev chain, as well as sus¬ 
ceptibility [65] and pair correlation [59] in long-range cou¬ 
pled superconducting fermions. However, these works 
have not provided a systematic comparison of multi¬ 
ple quantities between topological/trivial phases or upon 
topological transitions within a single frame (model). 

In this paper, we study various physical quantities 
of ID one-component fermions having proximity-induced 
pairing gap and interparticle short-range interaction, as 
a generalization of the Kitaev model. These quantities 



2 


are obtained from density-matrix-renormalization-group 
(DMRG) [79, 80] calculations and categorized into three 
groups: (i) topological properties, including ground-state 
degeneracy and entanglement spectrum; (ii) thermody¬ 
namic properties, including compressibility and suscep¬ 
tibility; and (hi) condensate properties, including pair- 
condensate fraction and Cooper-pair size. We take the 
topological region indicated by the first group as a refer¬ 
ence and investigate the behavior of the second and third 
ones in the parameter space. As a result, we shall find 
alternative signatures to identify the topological states 
and topological transitions. In addition, by tracking the 
multiple signatures, we show the topological phase dia¬ 
gram as a function of interaction. Finally, we study the 
behavior of several signatures against finite size, disor¬ 
der, as well as inhomogeneity, and discuss their stability 
under these conditions. 

The paper is organized as follows. In Sec. II, we in¬ 
troduce the model Hamiltonian and define the physical 
quantities of interest. Section III shows the setup of 
DMRG calculations. We present the results and discus¬ 
sions in Sec. IV. Section V is the conclusion. 

II. MODEL AND PHYSICAL QUANTITIES 

In this section, we present the model in considera¬ 
tion and physical quantities of interest. We consider 
one-component fermions in chain lattices where particles 
can scatter through the Gooper channel (or form Cooper 
pairs) due to a combined effect of external proximity- 
induced pairing and internal short-range interaction. If 
there is only the external effect, the Hamiltonian is of the 
Kitaev form [9], 

L-l 

Hk = Ho + J2 + H-C.) , (1) 

i=i 

where c| creates a fermion on site j, A' (taken real 
without the loss of generality) describes the proximity- 
induced pairing, L is the total number of lattice sites, 
and 

L-l L 

= (c]cj+i + H.c.) - (2) 

i=i i=i 

is the noninteracting Hamiltonian with number operator 
fij = c^jCj, nearest-neighbor tunneling t, and chemical 
potential If there is only the internal interaction effect, 
the Hamiltonian reads 

L-l 

= + (3) 

1=1 

where negative (positive) V represents attractive (repul¬ 
sive) interaction. In realistic cases of a nanowire [36-41] 
or an atomic chain [19-21] on a superconducting sub¬ 
strate, or a tube of quantum gases in higher-dimensional 
optical lattices [35], both effects can take place. We hence 


model the system with a phenomenological parameter 7 
to describe the relative strength of the two effects and 
write down our model Hamiltonian, 

FT = (1 - 7)iFK + 7 H 1 

L-l 

= Ho + '^ A -f H.c.^ -k Vfijfij+i, (4) 

i=i 

with independent variables A = (1 — 7 ) A' and V = jV, 
which can be self-consistently determined by microscopic 
degrees of freedom and/or realistic parameters of the sys¬ 
tem. (Below we aim to study physical characteristics of 
iF in a range of its parameter space rather than determine 
the parameters for a specific situation.) The Hamiltonian 
always conserves the even/odd parity of total number of 

particles N = (TV) = due to [FT, (—1)'^] = 0 and 

conserves TV (meaning [FT, TV] = 0) in the limit A —0. 

We will study the behavior of three groups of phys¬ 
ical quantities of interest in the topological and non- 
topological regions as well as upon the topological tran¬ 
sition. The first group, called topological quantities, in¬ 
cludes ground-state degeneracy and entanglement spec¬ 
trum degeneracy. The former one can be characterized 
by the energy difference 

6E = [Fieven “ Fiodd | (5) 

between the lowest eigen energies Fieven/odd of the 
even and odd blocks of the Hamiltonian [where 

((—I)^)even/odd = il]j respectively. The ground-state 
degeneracy 6E = 0 occurs as the manifestation of Ma- 
jorana fermions in the topological region and does not 
otherwise. 

The entanglement spectrum is a series of eigenvalues 
of reduced density matrix pR obtained by tracing out 
half spatial degrees of freedom of the ground-state wave 
function [V'g), 

PR = Trj<L /2 14) (41 > ( 6 ) 

where we consider L an even number and choose the 
ground state of the even-particle-number block when the 
even-odd degeneracy occurs {6E = 0). We obtain a set 
of 2N\ largest eigenvalues of pR, denoted by {Aj} and 
sorted as Xj > Xk if j < k, and compute the difference 
between paired elements defined by 

Nx 

(5A = ^ IA 2 J -1 - A 2 j|. (7) 

T=i 

Given sufficiently large Nx, the condition dA —0 guar¬ 
antees twofold degeneracy in the entanglement spectrum, 
which is contributed by a pair of Majorana fermions sep¬ 
arated in each entanglement subsystem and hence re¬ 
garded as a signature for their existence (see details in 
Refs. [81-83]). 

The second group involves thermodynamic quantities 
that can be expressed as the derivative of energy density 
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E = {H) jL with respect to the system parameters. We 
focus on the first derivative — which is equal to the 
particle density p = N/L from the Hellmann-Feynman 
theorem, and second derivatives — ■> 'which describe 

compressibility {X = Y = p.) or various susceptibility 
(X = Y = A or X = p,Y = A) of the system. A 
discontinuity or a peak of energy derivatives represents a 
drastic change in the behavior of the ground-state energy, 
which implies a cross between two lowest-energy states of 
the system. Such a cross is a necessary condition for a 
topological transition. 

The third group includes quantities describing conden¬ 
sate properties. Superconducting fermions can be con¬ 
sidered as a pair condensate in which a two-body or 
Cooper-pair state is macroscopically occupied [84, 85], 
analogous to the Bose-Einstein condensation in bosonic 
systems [86, 87]. The condensate properties are well char¬ 
acterized by the pair density matrix , whose element 
is defined in spatial coordinates as 


We have shown the model Hamiltonian and all of the 
physical quantities of interest. Among these quantities, 
twofold degeneracies of ground-state energy and entan¬ 
glement spectrum can be regarded as direct signatures for 
the topological states and Majorana fermions. They have 
also been proved for interacting systems [62, 82]. A large 
Cooper-pair size is directly related to the long-range en¬ 
tanglement for the presence of Majorana fermions, while 
the interaction effects on it are to be investigated. The 
energy derivatives such as compressibility and suscepti¬ 
bility can imply a cross of the two lowest-energy states, 
but its relation to the topological phase transition needs 
to be confirmed by the comparison with the direct sig¬ 
natures. Below we will take the ground-state energy 
and entanglement spectrum as a reference to pinpoint 
other signatures displayed in the thermodynamic quan¬ 
tities and/or condensate properties. In the next two sec¬ 
tions, we present the numerical setup, results, and dis¬ 
cussions. 


Pjfcp'fe'= (c]4cfe'C/)- (8) 

When the condensation occurs, pP'^*'' has an eigenvalue 
[~ 0(A)] largely compared to the others [~ 0(1)], 
representing the macroscopic occupation or the number 
of condensed pairs. This number defines the condensate 
fraction as 


P = 


xr" - 2 

N ’ 


(9) 


with an offset 2 of the noninteracting limit [60] [such 
that P = 0 for a free system described by iJg in Eq. (2)]. 
The eigenstate corresponding to Ag^^”^ represents the 
Cooper-pair wave function. The Cooper-pair size is de¬ 
fined as a root-mean-square distance between the two 
particles in a pair, 


(j - ky 




( 10 ) 


and indicates a length scale over which two fermions bind 
in real space. The topological state is a weak-pairing 
state [88] with rpair —t oo (or ^ L in a finite-size sys¬ 
tem), which originates the long-range entangled Majo¬ 
rana fermions on both ends. Instead, the trivial state has 
finite rpair (or ^ L). Therefore, a drastic change in rpair 
is expectable upon the topological transition. Note that 
the commonly used U{\) symmetry-breaking order [6] 
can distinguish the weak-pairing state (being constant, 
corresponding to an infinite Cooper-pair size) from the 
strong-pairing one (exponentially decaying in real space, 
with the decay rate determining a finite Cooper-pair size) 
for a number-nonconserving state, but fails to do so with 
a number-conserving one (e.g., A —^ 0 in our model). 
The definition of rpair in Eq. (10) works for both. (In fact, 
such rpair iias been applied to study topological proper¬ 
ties in a two-dimensional system [59].) 


III. NUMERICAL SETUP 

While exact solutions can be found in specific param¬ 
eter regions [89], the many-body ground state of the in¬ 
teracting Hamiltonian in Eq. (4) can, in general, be com¬ 
puted numerically. Our numerical results are obtained 
using the DMRG [79, 80] method, whose accuracy has 
been demonstrated in computing the ground-state prop¬ 
erties of a short-range coupled ID system. This method 
has been widely adopted to study topological proper¬ 
ties in spinless fermions [63, 90], spin-orbit coupled elec¬ 
trons [62, 66, 67], and ultracold atoms [91, 92]. In our 
work, we employ the DMRG method on systems up to 
L = 256. The Z 2 symmetry of parity conservation is con¬ 
sidered to reduce the computational cost. We keep up to 
m = 120 states and apply seven sweeps in the ground- 
state calculation. The number of states m kept in the 
DMRG method determines the size of the approximated 
Hamiltonian and hence the accuracy of calculations. The 
discarded weight in the eigenvalue of the reduced density 
matrix with m = 120 is of the order of 10“®, which guar¬ 
antees the convergence of the ground-state properties. 

The DMRG calculations are numerically efficient for 
obtaining the ground-state energy, entanglement spec¬ 
trum, and real-space two-point correlations. However, 
the bottleneck of our numerical study lies in the com¬ 
putation of the condensate properties. The construction 
of the pair density matrix requires the evaluation of all 
possible four-point correlations. Due to the L^-growing 
computational cost as the system size increases, we are 
bounded to L = 64. Another limit is the diagonalization 
of the pair density matrix, whose computational effort 
scales as L® and will eventually dominate in a large-size 
case. 


IV. RESULTS 

In this section, we present and discuss results that 
track multiple quantities for several cases of interest. 
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We take the tunneling strength t to be the energy units 
for convenience and set A = 0.2. First, we study the 
well-known Kitaev chain (where the internal interaction 
V = 0) and find the compressibility and susceptibility 
as good signatures for the topological transition. Sec¬ 
ond, we track the compressibility for general cases with 
the internal interaction. We will show that it remains 
a valid indicator because it indicates the same topolog¬ 
ical region as the entanglement spectrum does. Third, 
we investigate the condensate properties together with 
the compressibility as topological signatures for a rela¬ 
tively small system. We also study the characterization 
of topological phases with the use of number density and 
condensed-pair density. Finally, we present the effects of 
finite size, disorder, and inhomogeneity. 


A. The Kitaev model 

The Kitaev chain model is described by the Hamilto¬ 
nian of Eq. (1). It also represents the weakly interact¬ 
ing limit F —?► 0 of our model Hamiltonian of Eq. (4). 
The upper and lower topological transition points of an 
infinite Kitaev chain are known as /r = ±2 = [9], 

between (outside) which the system is in a topological 
(trivial) state. Eor a hnite chain, we numerically calcu¬ 
late the transition points. They show little change if the 
system size is large compared with the characteristic size 
of Majorana fermions. 

Eigure 1 shows four thermodynamic quantities as en¬ 
ergy derivatives: density p = —[Fig. 1(a)], compress¬ 
ibility [Fig. 1(b)], and susceptibilities — 

[Fig. 1(c)] as well as [Fig. 1(d)], vs the chemical 

potential p. We plot data for an inhnite chain (solid 
curve) and a hnite chain of L = 256 with open (red 
crosses), periodic (blue circles), and antiperiodic (green 
triangles) boundary conditions. We see that the density 
smoothly varies from vacuum (p = 0) to commensurate 
Filing (p = 1) as p increases. The trivial states have ei¬ 
ther a very low Filing (p ^ 0 in /r < p~) or an almost 
commensurate Filing (p ^ 1 in p > /i+). All of the sec¬ 
ond derivatives of the energy develop peaks at each tran¬ 
sition point. The compressibility peaks reFect the bulk 
gap closing at the transition point [78]. Across the peak, 
a sign change of the curve slope indicates a discontinuity 
in the third derivative of energy, which implies that the 
topological transition is third order. Such a topological 
transition type has also been found in a long-range cou¬ 
pled system [65]. The peak structure is independent of 
the chain size (see more analyses in Sec. IV D) as well 
as the boundary conditions and can thus be taken as a 
reliable signature for the topological transition. While 
the two susceptibilities both ought to work well, below 
we study the compressibility, i.e., a typical observable in 
experiments, as an indicator for cases with internal in¬ 
teraction. 



FIG. 1. (Color online) Noninteracting system energy’s (a) 
first derivative = P and various second derivatives (b) 

-0 = If’ ''S chemical 

potential p. Four cases of L = oo (solid curve) and L = 256 
with open (red crosses), periodic (blue circles), and antiperi¬ 
odic (green triangles) boundary conditions are presented. The 
superconducting gap is set as A = 0.2. 


B. Effects of internal interaction 

With the internal interaction turned on (V yf 0), we 
calculate the ground state properties as a function of V 
for a case of L = 256. In Fig. 2(a), we plot compressibility 
vs chemical potential p at various V = —0.4 (purple 
diamonds), —0.2 (blue circles), 0 (green triangles), 0.2 
(yellow squares), and 0.4 (red crosses). The curves main¬ 
tain the two-peak structure, while the right transition 
point p^ has a positive (negative) shift at repulsive (at¬ 
tractive) internal interaction, and the left one p~ shows 
little change with the interaction. We also confirm that 
the region sandwiched by the compressibility peaks coin¬ 
cides with the topological region identified by the ground- 
state degeneracy and entanglement spectrum. Therefore, 
compressibility can be regarded as a reliable signature for 
topological states of interacting systems. Figure 2(b) 
shows compressibility vs density p in the same conven¬ 
tion as Fig. 2(a). Similarly, the curves display peaks at 
the two transition points. We notice that the position 
of the left (right) peak is at a low (high) filling, p < 0.1 
(> 0.9), and is insensitive to interaction. In other words, 
the topological state survives in most intermediate-filling 
regions, which promises the Majorana fermions in a wide 
range of density-controllable systems such as ultracold 
atoms. 

Since the topological region shrinks as the interaction 
becomes more attractive, we turn to study the fate of the 
topological state in the strongly attractive region. (Note 
that the strongly repulsive region has been well studied 
in Ref. [63].) The top panel of Fig. 2(c) shows compress¬ 
ibility vs p at V = —1.5 (red crosses), —2 (green trian¬ 
gles), and —2.5 (blue circles), with the topological region 
shaded. We see that the two peaks at p^ merge into one 
at F = —2.5 and the topological region disappears. The 
disappearance can also be seen in the entanglement spec- 
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FIG. 2. (Color online) (a),(b) Compressibility ^ vs fj, and p, respectively, at various interactions V = —0.4 (purple diamonds), 
—0.2 (blue circles), 0 (green triangles), 0.2 (yellow squares), and 0.4 (red crosses) for an open chain of L = 256 and A = 0.2. (c) 
(top panel), p (middle panel), and entanglement spectrum (ES, bottom panel) vs p at relatively large attraction V = —1.5 
(red crosses), —2 (green triangles), and —2.5 (blue circles). The topological region for each V is shaded in the top and bottom 
panels for comparison. The dashed curve in the top panel shows how the two compressibility peaks would emerge at V = —2.5, 
although the actual value is undefined due to the discontinuity in p. (d) Interacting phase diagram showing the topological 
region in the V-p plane. The upper and lower boundaries are marked by red crosses and blue circles, respectively, while the 
dashed curve shows the Hartree mean-held results. 


trum in the bottom panel as the doubly degenerate region 
vanishes at V = —2.5. The density curves in the middle 
panel show that the topological transitions at V = —1.5 
and —2 still happen at either a low or high filling. When 
the topological region disappears at V = —2.5, the den¬ 
sity (= curve exhibits a discontinuity and the sys¬ 

tem undergoes a first-order transition between a trivial 
low-filling state and a trivial high-filling one. 

The phase boundary shift can be explained by an ef¬ 
fective chemical-potential shift due to the interaction. 
Considering a Hartree approximation rij+irij —>■ pnj + 
prij^i — , one can turn the Hamiltonian of Eq. (4) into 

the original Kitaev form with an effective chemical po¬ 
tential /Xeff = p — ‘2Vp, where p = p{p,V). The upper 
and lower transition points are thus given by pes = ±2, 
respectively, which leads to 

p^ =±2 + 2Vp{pf,V). (11) 

Because the density at the lower (upper) transition point 
is low (high) and insensitive to interaction, we approx¬ 
imately insert p{p~,V) = 0 and p{p'^,V) = 1 into 
Eq. (11) and obtain = 2 -|- 2E and p~ = —2. Such 
relations tell that the lower boundary barely depends on 
V and the upper boundary linearly shifts with V. Figure 
2(c) shows the topological region in the V-p plane. The 
upper and lower phase boundaries from this approxima¬ 
tion (dashed curves) well match those from the numerical 
calculations (red crosses and blue circles, respectively). 

C. Condensate properties 

In this section, we study pair-condensate properties of 
the system. As mentioned in Sec. HI, the calculation of 
the pair density matrix can be time consuming. Such 


a constraint directs our focus on a relatively small sys¬ 
tem (L = 32) rather than a large one. The top panel 
of Fig. 3 shows the compressibility curves as a reference 
for the topological region at various interactions [con¬ 
ventions are the same as in Fig. 2(a)]. The middle and 
bottom panels show the condensate fraction P defined 
in Eq. (9) and Cooper-pair size rpair defined in Eq. (10), 



FIG. 3. (Color online) Compressibility (top panel), con¬ 
densate fraction P (middle panel), and Cooper-pair size Tpair 
(bottom panel) vs p at various interactions for an open chain 
of L = 32 and A = 0.2. Conventions are the same as in 
Fig. 2(a). 
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respectively. Different from the compressibility and en¬ 
tanglement spectrum, we first see that the curves of con¬ 
densate properties lack symmetry with respect to the 
half-filling point. At the upper phase boundary (right 
peak of the compressibility curve) both P and rpair de¬ 
velop a kink, indicating a sudden change in the con¬ 
densate properties upon the topological transition. The 
sharp decrease of rpair is consistent with the transition 
from weak-pairing (topological) to strong-pairing (triv¬ 
ial) states [6, 88]. At the lower boundary (left peak of 
the compressibility curve), rpair shows a peak, while P 
changes the trend but does not show a clear signature 
due to the finite-size effect (the filling is so low such that 
the total number of particles is smaller than two). We in¬ 
crease the system size and find that the turning point of 
P approaches the phase boundary. Therefore, the behav¬ 
ior of P and rpair can be an indicator for the topological 
transition. 

The topological transition accompanied by a signifi¬ 
cant change in the Cooper-pair size raises an interesting 
question of whether the change in the pair condensate 
can be purely qualitative or must be both qualitative and 
quantitative. We try to answer this question by exam¬ 
ining two essential quantities of a pair condensate, i.e., 
the total number of particles N and the total number 
of condensed pairs around the topological phase 

boundary. The answer is the former if the system can go 
across the phase boundary and keep both N and 
unchanged (only rpair changes). Our model is suited to 
explore this question because it is defined in a three- 
dimensional parameter space {A,V,^) such that a path 
along which two functions N{A,V,fi) and Aq‘*'”^(A, V, /i) 
remain constant becomes mathematically possible. We 
explore large enough regions around upper and lower 
transition points, {A,V,fi) = (0.2,0, ±2), respectively, 
by varying all three parameters for a L = 64 case. The 
results show that one of N and can remain constant 
upon the topological transition, but not both. In other 
words, the topological transition or the sudden change of 
the Cooper-pair size in our model can be regarded as a 
result of a quantitative change in N or Moreover, 

we do not find that a topological state and a trivial state 
have the same N and Ag'^”. Our results may have impli¬ 
cations for the characterization of topological states by 
(iV/L, AgP^VL). 

We further study this possibility by computing the con¬ 
densate properties for hundreds of points in the range of 
0.15 < A < 0.45, -0.45 < D < 0.4, and 1.4 < /r < 3.4, 
for L = 32. Figure 4 shows the data points representing 
three different ranges of the Cooper-pair size around the 
upper topological transition region in the plane of N/L 
and Ag“'^/L. The system can be considered as a topo¬ 
logical state if Tpair > 4. We see that along a vertical 
(horizontal) path in Fig. 4, the system can undergo a 
topological transition at fixed N (Ag'^'’^). On the other 
hand, the concurrence of a decrease in rpair and an in¬ 
crease in either N or Ag*^'*^ confirms our conjecture, i.e., 
an inevitable quantitative change upon the topological 



N/L 

FIG. 4. (Color online) Data points showing three different 
ranges of the Cooper-pair size, i.e., rpair < 3 (blue circles), 
3 < rpair < 4 (red crosses), and rpair > 4 (green triangles), in 
the plane of N/L and for L = 32. 

transition. In other words, rpair is a function of N and 
Ag*^**^, and these two parameters can hence be used to 
characterize the topological phase diagram. We finally 
comment that a more solid confirmation lies in the con¬ 
vergence of Fig. 4 in the thermodynamic limit; we leave 
the test of this for future study once more powerful com¬ 
putational tools are available. 

D. Effects of finite size, disorder, and 
inhomogeneity 

In the last section, we explore the stability of the signa¬ 
tures against three realistic effects in experiments—finite 
size, disorder, and inhomogeneity. Figure 5 shows the 
compressibility (top panel), energy difference SE (mid¬ 
dle panel) defined in Eq. (5) , and entanglement spectrum 
difference 5X (bottom panel) defined in Eq. (7) vs fj, for 
L = 32 (red crosses), 48 (blue circles), and 64 (green tri¬ 
angles). We see that although all three quantities indi¬ 
cate the transition points, the compressibility curve is in¬ 
sensitive to the system size, while SE and 6X can increase 
by several orders as the system size is halved. Therefore, 
the compressibility is a stable indicator against the finite- 
size effect. The insensitivity also provides a more numer¬ 
ically efficient way to predict the topological transition 
of a large-size system by calculating the compressibility 
peak of a relatively small one. 

For the disorder effect, we consider a Hamiltonian 

L 

Hd = H + Y,Sjnj, ( 12 ) 

1=1 

with a set of random local potential shifts {(5^ } that obey 
the normal distribution and have zero average. For the 
inhomogeneous effect, we consider an external harmonic 
trap with curvature K turned on. 
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FIG. 5. (Color online) Compressibility (top panel), en¬ 
ergy gap between the ground state and first excited state SE 
(middle panel), and difference between even and odd sectors 
of entanglement spectrum dA (bottom panel) vs p around 
the lower (left column) and upper (right column) transition 
points at various open-chain sizes: L — 32 (red crosses), 48 
(blue circles), and 64 (green triangles), respectively. Data are 
for A = 1/ = 0.2. 



FIG. 6. (Color online) Comparison between a disorder system 
[Hd of Eq. (12), red crosses], a trap system [Ht of Eq. (13), 
blue circles], and a reference without these effects [H of 
Eq. (4), green triangles]. Conventions are the same as in 
Fig. 5, except the system size here is L = 128. 


In Fig. 6, we compare a disorder case (red crosses) with 
the variance of {6j} equal to 0.1 and a trap case (blue cir¬ 
cles) with K = 0.8/(L—1)^ with the original Hamiltonian 
H ioi L = 128 in the same convention of Fig. 5. We see 


that neither of the two effects can alter the signatures for 
the topological transition. Such results are anticipated 
because the topological states and Majorana fermions are 
symmetry protected. Perturbations can not destroy the 
topological order as long as they are not strong enough 
to cause the bulk-gap crossing. (Regimes of strong dis¬ 
order or deep traps are beyond the scope of this study. 
One could refer to previous works in Refs. [34, 93-96].) 

V. CONCLUSION 

In conclusion, multiple physical quantities have been 
analyzed for one-component fermions with proximity- 
induced superconducting gap and interparticle interac¬ 
tion in ID lattices, which can be a topological supercon¬ 
ductor hosting Majorana fermions. In addition to the 
double degeneracy of ground-state energy and entangle¬ 
ment spectrum, we have found that the topological tran¬ 
sition can also be revealed by peaks of compressibility and 
susceptibility curves, as well as a sudden change of trend 
in the condensate fraction and Cooper-pair size. Among 
them, the compressibility peak is particularly useful for 
its stability against the finite-size effect and being ob¬ 
servable in experiments. The Cooper-pair size directly 
shows the topological transition between strong-pairing 
and weak-pairing states. By tracking these signatures, we 
have found that the topological transition is third order. 
As the interaction becomes more attractive, the topolog¬ 
ical state finally disappears and the system undergoes a 
first-order transition between low-filling and high-filling 
trivial states. We have also explored the possibility to 
characterize the topological phase using density of parti¬ 
cles and that of condensed pairs. One future direction is 
the extension of this study to other interacting platforms 
in which various tunneling or pairing channels [60, 97, 98] 
need to be considered and an alternative treatment suited 
for the continuous space [99, 100] may apply. In addition, 
our results may find applications on spin systems that 
are associated with our Hamiltonian of Eq. (4), such as 
an Ising chain with a transverse field {V —>■ 0, e.g., see 
Ref. [101]), the XXZ model (A —>■ 0) [102], the Bax¬ 
ter XYZ model —?► 0) [103], and the two-dimensional 
classical Ising model (U —?► 0, t = A) [104]. 
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